Quantum Zeno and anti-Zeno effects in an Unstable System with Two Bound State 
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We analyze the experimental observations reported by Fischer et al [in Phys. Rev. Lett. 87, 
040402 (2001)] by considering a system of coupled unstable bound quantum states \A) and \B). The 
state \B) is coupled to a set of continuum states \CQ(u))). We investigate the time evolution of |j4) 
when it decays into |CO(oj)) via \B), and find that frequent measurements on \ A) leads to both the 
quantum Zeno effect and the anti-Zeno effects depending on the frequency of measurements. We 
show that it is the presence of \B) which allows for the anti-Zeno effect. 



PACS numbers: 03.65.Xp, 03.67.Lx 

I. INTRODUCTION 

The quantum Zeno effect, first predicted by Misra and 
Sudarshan [H, , is the hindrance of the time evolution 
of a quantum state when frequent measurements are per- 
formed on it. In the limit of continuous measurement 
the time evolution of the state, in principle, completely 
stops. The seminal paper by Misra and Sudarshan proves 
the existence of an operator corresponding to continuous 
measurement belonging to the Hilbert space of a generic 
quantum system. More recently, several authors have 
suggested that the opposite of quantum Zeno effect may 
also be true [1,0, [f|. That is, frequent measurements can 
be used to accelerate the decay of an unstable state. This 
effect is known as the anti-Zeno effect or the inverse Zeno 
effect. The original formulation of the quantum Zeno 
effect treated the measurement process as an idealized 
von-Neumann type; that is an instantaneous event that 
induces discontinuous changes in the measured system. 
The anti-Zeno effect was first identified as a possibility 
when measurement processes that take a finite amount 
of time were considered. This led to the suggestion by 
several authors that the anti-Zeno effect should be ob- 
served more often in physical systems than the quantum 
Zeno effect. 

Experimental evidence supporting the quantum Zeno 
effect in particle physics experiments was first pointed 
out by Valanju et al @, 0|- Direct experimental observa- 
tion of the quantum Zeno effect was obtained by Itano 
et al [8] in a three-level oscillating system. Recently, 
in a set of experiments Fischer, Gutierrez-Medina, and 
Raizen observed, for the first time, both the quantum 
Zeno effect and the anti-Zeno effect in an unstable quan- 
tum mechanical system In this Letter we present a 
simple model that reproduces all of the important results 
of this experiment. 

This Letter is organized as follows: We briefly describe 
the experiment by Fischer et al in section [TTJ In section 
IIII1 we present a simplified model of the system studied 
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experimentally in Q . The model is exactly solvable and 
the time dependence of the survival probability of the 
initial unstable state can be analytically calculated. In 
section ITVl we show that the solutions reproduce all the 
important features of the experimental system. On the 
basis of our model we argue that the anti-Zeno effect 
is observed because of the presence of more than one 
unstable bound states in the system. Our conclusions 
are in section [Vj 



II. DESCRIPTION OF THE EXPERIMENT 

In the experiment by Fischer et al 9], sodium atoms 
were placed in a classical magneto-optical trap that could 
be moved in space. The motional states of the atoms in 
the trap were studied. Initially, the atoms were placed in 
the "ground" state of the moving trap so that they re- 
mained inside the trap. The stable bound states occupied 
by the atoms in the trap are made unstable by acceler- 
ating the atoms along with the trap at different rates. 
By tuning the acceleration appropriate conditions are 
created for the atoms to quantum mechanically tunnel 
through the barrier into the continuum of available free- 
particle states. The number of atoms that tunneled out 
of the bound state inside the trap as a function of time 
was estimated at the end of the experiment by recording 
the spatial distribution of the atoms. Since the trap was 
accelerated throughout the experiment, atoms that spent 
more time in the and trap had higher velocities and they 
moved farther in unit time. So by taking a snap-shot of 
the spatial distribution of all the atoms at the end of the 
experiment the time at which each one of them tunneled 
out of the trap can be estimated. 

To obtain the Zeno and the anti-Zeno effects, the tun- 
neling of the atoms out of the trap has to be interrupted 
by a measurement that estimates the number of atoms 
still inside the trap. Such a measurement was imple- 
mented in the experiment by abruptly changing the ac- 
celeration of the trap so that tunneling from the ground 
state is temporarily halted. These interruption periods 
were long enough (AOfis) to separate out the atoms that 
tunnel out before and after each interruption into resolv- 
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FIG. 1: On the left, the lower line is the "unmeasured" decay 
curve corresponding to the case where the trapped atoms are 
accelerated with no interruptions so that tunneling out of the 
trap is always present. The upper line corresponds to the 
case where the tunneling is interrupted every 1 /is leading 
to the quantum Zeno effect. On the right, the upper line is 
the "unmeasured" decay and the lower line corresponds to 
interruptions every 5 /xs leading to the anti-Zeno effect in the 
experiment by Fischer et al. 

able groups. By measuring the number of atoms in each 
group and knowing the total number of atoms that were 
initially in the trap, the number of atoms in the trap 
at the beginning of each interruption period was esti- 
mated. Using this data, the time dependence of the sur- 
vival probability of the bound motional states of an atom 
was reconstructed. Their observations are summarized in 
Fig. [TJ 

The frequency with which interruption periods were 
applied determined whether the Zeno or anti-Zeno ef- 
fects were obtained. Repeatedly interrupting the system 
once every micro-second led to the quantum Zeno effect. 
With an interruption rate of 5/is the anti-Zeno effect was 
obtained. 

In FhjJTJ the zero slope for the survival probability at 
t = is expected for the time evolution of a generic 
unstable quantum state on the basis of the analyticity of 
the survival probability and its time reversal symmetry 
(Tfj| . This non-exponential behavior at short times leads 
to the quantum Zeno effect. 

The shape of survival probability as a function of time 
of the "unmeasured system" has an inflection point at 
t 7fis. (see Fig. [T]). We show that it is the presence 
of a second unstable bound state that is responsible for 
this inflection point. The anti-Zeno effect is obtained by 
taking advantage of that inflection point. 

In an earlier experiment that led to this experiment, 
Bharucha et al. observed tunneling of sodium atoms from 
an accelerated trap [ll[ . Our analysis of the experiment 
in Q is motivated by the following passage from [TT| : 

When the standing wave is accelerated, 
the wave number changes in time and the 
atoms undergo Bloch oscillations across the 
first Brillouin zone. As the atoms approach 
the band gap, they can make Landau-Zener 
transitions to the next band. Once the atoms 
are in the second band, they rapidly undergo 



transitions to the higher bands and are effec- 
tively free particles. 

The last sentence above suggests that higher energy 
bound states are present in their system. An atom in the 
ground state might have to go through these intermediate 
bound states before it can tunnel out to the continuum 
of available free-particle states. We investigate the effect 
of these intermediate states on survival probability of the 
ground state. 

III. MODEL 

We consider an interacting field theory of four fields 
labeled A, B, C, and O, with the following commutation 
relations: 

[a,at] = [6,6t] = [ c , c t] = i, 

[6(u),tf(u')j\=6(u>-u/). 

All other commutators being zero, (a) etc. repre- 
sent the creation (annihilation) operators corresponding 
to the four fields. Only 9 is labeled by a continuous 
index u, while other fields are assumed to only have dis- 
crete modes. The allowed processes in the model are 

A ^ B and B <=> CO. 

The Hamiltonian for the model with these allowed pro- 
cesses can be written down as 

H = H + V (1) 

where, 

H Q = E A a i a + E B tfb+ duj u <9 f {uj)6{uj) (2) 
Jo 

and 

V = at b + fl* a 

+ du [f(u) c 6(u) + f{uo)* J6\oj) b] . (3) 
Jo 

The two discrete energy levels are denoted by Ea and 
Eb- The Hamiltonian in Eq. (fT]) is obtained by modifying 
the Hamiltonian for the Friedrichs-Lee model [l^, EH . It 
is instructive to look at the spectrum of H and H in the 
complex plane at this point. We call the eigenstates of H 
physical states, while the eigenstates of Ho are referred 
to as bare states. We choose the zero point of energy 
so that the continuum eigenstates of Hq have positive 
energy (0 < lo < oo). If Ea and Eb are negative and 
if the shift in these energies due to the perturbation V 
is small then the spectrum of the physical Hamiltonian, 
will contain two stable bound states with negative 
energies A a and A^ in addition to continuum of states 
with positive energies. This is illustrated in Fig. O 
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FIG. 2: Ea,b are the energies of bare bound states. We choose 
Ea and Eb to be real and negative. The full Hamiltonian, H, 
then has two real negative eigenvalues indicated by Aa and 
Ab ■ H also has a continuum of eigenstates along the positive 
real- axis. 
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FIG. 3: The discrete eigenvalues of Ho, Ea.b are chosen to be 
real and positive. The full Hamiltonian has no real negative 
eigenvalues corresponding to bound states. The resolvent of 
H has two complex poles. The two poles represent unstable 
states and they could have either positive or negative imagi- 
nary parts. These four possibilities are indicated by Ai,2,3,4- 
The sign of the imaginary part is fixed by the boundary con- 
ditions. A2 and A4 correspond to unstable states that decay 
in time. The continuum of physical states lies along positive 
real- axis. 



We are interested in studying the temporal evolution of 
an unstable system. So we choose Ea and Eb to be posi- 
tive so that they lie embedded in the physical continuum 
spectrum. The spectrum of the physical Hamiltonian 
will no longer include bound states. The eigenstates of 
H belonging to the continuum, corresponding to eigen- 
values < A < 00, will form a complete set of states. To 
see what happened to the two bound states of the bare 
Hamiltonian Hq when the perturbation V is introduced, 
it is instructive to look at the resolvent of the physical 
Hamiltonian, (E — H)^ 1 . The resolvent, in this case, 
will have two complex poles indicating two transient or 
unstable states. The location of these poles is fixed by 
choosing appropriate boundary conditions to make sure 
that the unstable states decay (rather than grow) in time. 
This is illustrated in Fig. [3J 



A. The Continuum States 

We are interested in the time evolution of the eigen- 
states of Hq, namely the two bound bare states \A) and 
\B), and the continuum states \CQ(lu)). The state \A) 
in our model corresponds to the unstable bound state 
occupied by the atoms inside the trap in the experiment 
by Fischer et al The states \CQ{uj)) represent the con- 
tinuum outside the trap into which the bound state can 
decay. 

We have introduced an additional unstable bound state 
\B) which represents a second bound motional state of 
the trap. The state \ A) is directly coupled to only \B) and 
the decay of \A) into |CG) is mediated by the new state 
\B). We will show that the presence of the additional 
bound state can explain several of the key features of the 
experiment by Fischer et al 

We begin by writing the full Hamiltonian, H , in matrix 
form using eigenstates of Hq as basis [14| . 



H 



n* 

Eb /*(«/) 

f(uj) U>S(U) — Ld') , 



(4) 



Let \ip\) represent an eigenstate of H with eigenvalue A, 
satisfying the eigenvalue equation 

#V>A = AV-A- (5) 
We express ip\ also in terms the eigenstates of H , 



( Wa> 
Wa> 

V(ce(w)|^) 



(6) 



H can have three possible classes of eigenstates; a maxi- 
mum of two bound states and a set of continuum eigen- 
states. We first look at the continuum eigenstates of H 
followed by the remaining physical bound states in next 
section. 

Using the Eq. we get a system of three coupled 
integral equations: 



(A ? E A ), 



X-E a ' 

X-E B 



A-, 



(7) 

(A 7^ Eb), (8) 
(9) 



The continuum wave function in Eq. §§§ is singular in the 
energy (momentum) space, therefore it extends to infinity 
in the configuration space. We choose the "in" solution 
by choosing the sign of the imaginary part. Eqs. ([7][9]) 
can be solved simultaneously to get 
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, /(A) 

\ /3+(A) X-uj+ie 



/(A) n* 

/3+(A) \-E A 



(10) 
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where 

0(z) = z-E B 



O 2 



(11) 



(12) 



z-E A 

is a real analytic function and 

/3 ± (A) = /3(A ± ie) and /3 ± (A)* = /3 T (A). 

B. The Physical Bound States 



Now we look at the case when A ^ u> for all A. Since 
< lo < oo, the physical states in this case are restricted 
to eigenvalues on the negative real-axis. The delta func- 
tion is no longer present in Eq. @, meaning that the 
physical states do not extend to infinity. Such solutions 
of the eigenvalue problem correspond to the stationary 
eigenstates of H . Eq. ([9]) now becomes 



/M B 

A^/ A ' 



Solving for /i^, we get 



/3(AK = 0. 



(13) 



(14) 



We choose /if ^ 0, to obtain a non-trivial solution, then 
/3(A) = 0, which leads to the following expression, with 
zeroes at A = A,*. 



A 2 - A [E A + E E 



l/HI 2 



-n 2 + e a e 



A — ui 

x l/MI 2 

A — oj 



du> 



(15) 



The wave functions belonging to this class of solutions 
are just complex numbers that are fixed by the normal- 
ization condition. The solutions are the following: 



^A j 

where 




( ./Wi 



1 

i jM 



\ 



(16) 



dft_ 



(17) 



A=A, 



= 1 



{\-E A y 7o (A-w) 



: duj 



We can find Aj by numerically solving Eq. (TTS"]) . For a 
weak potential V, we can roughly approximate 



A 7 - ss E A b ± i Im 



l/MI 2 



A — lu + ie 



\=E A , 



(18) 



In general, Eq. (| 1 5|) has four possible roots, given by 
Eq. (TT8"]) . Only the bound states with real eigenvalues 
Aj are physically relevant. Once again this due to the 
fact that the physical spectrum is confined to the real- 
axis. The imaginary part of the integrals in Eq. (fl8|) 
vanishes everywhere except near luq. Hence for E A and 
Eb with sufficiently negative energies there are two phys- 
ical bound states with real eigenvalues. Once again this 
will not be the case of interest since no decay takes place 
here. For sake of completeness it should be noticed that 
the physical bound states will be present wherever the 
imaginary part of the integral in Eq. (|18[) vanishes. 



C. Survival Amplitude 

We are now in a position to calculate the survival am- 
plitude of the bare state as a function of time. If |^4) is 
occupied at t=0, then its survival amplitude is obtained 
by computing the matrix element £sf A (t) = (A\e~ lHt \A). 
We can compute this matrix element by inserting a com- 
plete set of physical states, i.e. 



IV'a) (V>a| d\ + |^A A ) (lpA A \ + \ll>A B ) (V'AbI = 1- 



Here we are considering the possibility that stable bound 
states \tpA A B ) w hh real energies K a _b may exist. Let's 
continuing with the calculation of the survival amplitude. 



s* A {t) = / {A\e~ Mt \^ x ) {ip x \A) d\ 



{A\e~ lHt 
+ (A\e 



\i>A A )&A A \A) 

- im \i>A B ) &A B \A) 



-ikt 



+e 



{A\i> x ) (ip x \A) dX 



-iA A t 



(4^a*)<Va a |A) 



+e-^ (A\ip AB ) (1>a b \A) 



lXt \vt\ 2 dX- 



e- iAAt \nt 



-e-^*|< 



(19) 



We have used the completeness and orhtonormality of the 
physical states in the above equation. Orthonormality of 
these states is shown in Appendic A, while completeness 
can be demonstrated by using the techniques described 
in Appendix A and in reference [l4l |. 

The survival probability of \A) is simply the square of 
the amplitude, 



P A (t) = KaWI' 



Ma 



dX 



-iA B t I n A 
MAi 



.(20) 
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Similarly, if state \B) is occupied at t=0, then the survival 
probability of \B) will take the form 



P B {t) = \^B{t)f 



-iXt I ,.B 



n"\ dX + e 



-iA A t\ ,.B |2 



(21) 



IV. 



NUMERICAL CALCULATIONS AND 
RESULTS 
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In this section, we discuss the numerical calculations 
of the survival probabilities Pa(*) and Ps(t). We choose 
the form factor to be 



(OJ - UJ ) 2 + /i 2 



(22) 



The factor of y/uj in the numerator is a phase-space con- 
tribution. The rest of f(u>) is just the Lorentzian line 
shape. The function f(u>) peaks near ujq and its width is 
controlled by /z, and its strength by a. 

We require the strength of the perturbation to be weak 
relative to the eigenvalues of Hq. In other words, we 
don't want the original system to change drastically. So 
we choose small values for the parameters fi, a and £1 
and set fi = 0.30, a = 0.11, and Q = 0.04. We also 
want the bare bound states to become unstable in the 
presence of the perturbation V. This can be achieved 
by setting their energies to lie in the physical continuum 
sufficiently above the threshold. The physical continuum 
ranges from zero to infinity, and so we choose numerical 
values E A = 2.00 and E B = 2.10. 

For the choice of Ea and Eb here, Eq. (fTS]) yields 
complex eigenvalues for the physical bound states. The 
condition on the physical bound states' stability requires 
that they have real-negative eigenvalues and so here they 
are unstable and show up only as a spectral density. Since 
the physical spectrum is confined to the real-axis and 
since there are no stable physical bound states, the last 
two terms in the Eqs. ([20]) and (f2Tj) are zero. The physical 
continuum states \ip\) form a complete set of states by 
themselves. The equations for survival probability then 
reduce to 



and 



P A {t) 



P B {t) 



-iXt 



Ht\ dX 



-iXt I ,.B 



Ma dX 



(23) 



(24) 



A. "Unmeasured" Evolution 

First, we study the survival probability as a function 
of time of \A) (Pa)- We then compare it to the survival 



FIG. 4: Survival probability of \A) given by Eq. fl23J|. The 
dashed line shows the "unmeasured" evolution of \A). The 
solid line shows the effect of repeated measurements made 
at high frequencies leading to the Zeno effect. The dotted 
line show the effect of repeated measurements made at lower 
frequencies leading to the anti-Zeno effect. The parameter 
values used are E A = 2.00, E B = 2.10, u) = 2.10, fi = 0.30, 
a = 0.11, and Q = 0.04. 
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FIG. 5: Survival probability of \B) given by Eq. (f24)l with 
(fi = 0). The dashed line shows the "unmeasured" evolution 
of \B). The solid line shows the quantum Zeno effect appear- 
ing due to repeated measurements. The parameter values 
that were used are Eb = 2.10, coo ~ 2.10, /i = 0.30, a — 0.11, 
and f2 = 0. 



probability of \B) (Pb) with Q — * 0. In the second case 
the first bound state |^4) is completely cutoff from the rest 
of the system as seen from Eq. @. Then we are deal- 
ing with a system with only one bound state coupled to 
the continuum. The differences between these two cases 
show precisely the effects that the second bound state 
has on the survival probability of \A). Note that the sec- 
ond case is nothing more than the well known Fredrichs- 
Lee Model. The form factor f(u>) is peaked around the 
Eb and so we have uio = 2.10. The "unmeasured" sur- 
vival probability of the states \A) and \B) are the dashed 
curves in Figs. 0] and O respectively. Notice that in both 
cases for long times the decay is exponential. 
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B. Effective Evolution: Zeno and anti-Zeno 



Interrupting the system by changing the acceleration 
is a necessary step in the experiment by Fischer et al to 
create the quantum Zeno effect or the anti-Zeno effect. 
The interruption resets the system and the tunneling pro- 
cess has to restart after the interruption period. In our 
model the shifting of is analogous to changing the ac- 
celeration in the experiment. The time evolution of the 
system (when SI = 0.04) can be hindered by shifting $1 
to a value much smaller than the difference between Ea 
and Eb- This effectively cuts off the oscillations between 
\A) and \B). During the interruptions, \B) is still con- 
nected to the continuum, just as in the experiment. The 
population of \A) remains constant during the interrup- 
tion periods. The length of the interruption period is 
important for resolving the different momentum states 
in the experiment by Fisher et al as shown in the last 
figure in @ . We do not have this constraint in our model 
and our measurements can treated as instantaneous (von- 
Neumann type). 

The experimental results presented in Q show the ef- 
fective evolution of the atoms in the trap with the in- 
terruption periods removed. We also look at the effec- 
tive evolution, with the interruption periods filtered out. 
Numerically, the procedure for obtaining the time evolu- 
tion with interruptions is simple. We start by fixing the 
interval of time, r, between the measurement induced 
interruptions. We start with a bare state with unit am- 
plitude and compute its survival porbability till time r. 
At this point the measurement is assumed to reset the 
system. The initial bare state wave function had only 
one non-zero component when expressed in the basis of 
bare states. Time evolution of this state under the full 
Hamiltonian makes all three components non-zero in gen- 
eral. Resetting the system correponds to setting the two 
new components that appeared as a result of the evolu- 
tion back to zero. This new (un-normalized) state is the 
starting point for further evolution until the next inter- 
ruption. This process is repeated several times to obtain 
the graph of the survival probability of the initial unsta- 
ble state when it is subject to frequent interruptions. 

The solid line and the dotted line in Fig. 2] shows how 
the effective time evolution of \A) can be hindered or 
accelerated by repeatedly interrupting the system. The 
second inflection point in the unmeasured evolution of 
the state \A) that allows us to obtain the anti-Zeno ef- 
fect is present due to the second unstable state \B). In 
the second scenario when we have set = there is only 
one unstable state in the system. From the graph of the 
survival probabilty of the unstable state without inter- 
ruptions shown by the dashed line in Fig. O we see that 
there is no way we can choose an interruption frequency 
that will lead to the anti-Zeno effect. On the other hand, 
interruptions made very frequently can lead to the quan- 
tum Zeno effect as shown by the solid line in Fig. 



V. CONCLUSION 

We present a simple theoretical model of the experi- 
ment done by Fischer et al In our model state \A) rep- 
resents the motional ground state in their experiment. 
\B) represents a higher energy motional bound state, 
while \C'Q(uj)) represents the continuum of available free- 
particle states in their experiment. The presence of the 
intermediate states is clearly stated in reference . We 
study the survival probability of the ground state as it 
decays into the set of continuum states via an intermedi- 
ate bound state. Our results are in agreement with the 
results of the experiment. 

We have shown how repeated interruptions of the time 
evolution of an unstable state can lead to the Zeno effect 
in our model. The anti-Zeno effect is also obtained in a 
similar fashion but we find that the system must have 
additional peculiarities for obtaining this effect. In our 
model there is an additional unstable state that mediates 
the decay of the original unstable state. The presence of 
this new state is shown to produce an inflection point 
in the graph of the survival probability of the unstable 
state that is of interest. This inflection point defines a 
time scale at which repeated measurements may be done 
on the system to effectively speed up its decay and thus 
obtain the anti-Zeno effect. We also point out that the 
second bound state is not needed to obtain the Zeno ef- 
fect. A generic decaying system can exhibit the Zeno 
effect if its time evolution is very frequently interrupted 
by measurements. We note here that a similar analysis 
has been done for oscillating systems by Panov llal a nd 
on the Friedrichs Lee model by Antoniou et al in [161 ] . 
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APPENDIX A: ORTHONORMALITY 

Here we want to show that the continuum states are 
orthonormal. That is 

I < \ 

= /if (^f,0AH) 

= MfMA+Mfuf + ^»0a(w) (Al) 
= S(X-r,). 
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We start by looking at the last term in Eq. (IA1[) . 

, /*fo)/(A) 
/3-(?7)/3 + (A) 

!« 

(A + ie — to) (77 — ie — w) 



/o 



+<5(A - 17) 



(A2) 



We can break up the integral in two integrals using par- 
tial fractions: 



1 



1 



(A + ie — lo)(t] — ie — uj) rj — X — lie 
1 1 



(A3) 



\ + ie — lu 77 — ie — ui 

Using Eq. (|A3[) . we re- write the integral term in Eq. (|A2 

as, 



l/HP 



1 



(A + ie — uj) (77 — ie — uj) 77 — A — lie 

-r*> |/( - )|2 i.(A4) 



/ A + ie — uj Jo 77 — ie — u> y 
Using Eqs. (fTTj) and (fT2)) , we can re- write Eq. (|A4|) in the 



following manner. 
1 



7/ — A — 2ie 



X + ie- j3+{X)- E B - 



-77 + ie + (3 (77) + E B + 



A - E A 



V-E A 
which then reduces to 

l^l 2 P~(V) /? + (A) 



-1- 



(A - E A )(rj - E A ) X + ie-r/ rj — ie — X' 



(A5) 



(A6) 



The last two terms in Eq. (|A6p cancel the first two terms 
in Eq. (|A2j) . and it reduces to 



r(r?)/(A) 



p-(r,)/3+(X)(r,-E A )(X-E A ) 



+ S(X-r,). (A7) 



Notice, the first two terms in the last equation are the 
negatives of y£ fxf + ^ ftf, and so the final result is 



r v <P\ = 6(X-r ) ). 



(A8) 



We have now shown that the set of continuum states is 
orthonormal. Using similar techniques, we can also show 
that the physical eigenstates form a complete set. 
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